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Abstract 

The coupled dynamics of the scissors mode and the isovector giant quadrupole 
resonance is studied in a model with separable quadrupole-quadrupole residual inter- 
actions. The method of Wigner function moments is applied to derive the dynamical 
equations for angular momentum and quadrupole moment. Analytical expressions 
for energies, B(M1)- and B(E2)-values, sum rules and flow-patterns of both modes 
are found for arbitrary values of the deformation parameter. Some predictions for 
the case of superdeformation are given. The subtle nature of the phenomenon and 
its peculiarities are clarified. 



1 Introduction 



The low-energy orbital magnetic-dipole excitations of deformed nuclei, commonly called 
the scissors mode, were predicted about 25 years ago [1, 2]. The prediction was inspired 
by the geometric picture that the (prolate) deformed neutron and proton distributions 
counter rotate (like scissors) within a small opening angle 5(j) in an oscillatory way around 
their common axis perpendicular to the long symmetry axis of the nucleus. Only a few 
years later the idea was confirmed experimentally with its detection in ^^^Gd [3]. At 
present, in addition to the rare earth nuclei, this excitation is also known in the actinides 
and in light nuclei. A complete review of the experimental situation can be found in 
[4]. The discovery of the scissors mode has initiated a cascade of theoretical studies. 
An excellent review of the present situation in this field is the one by D. Zawischa [5] 
(see also [6]). Very briefly the situation can be described in the following way. All 
microscopic calculations with effective forces reproduce experimental data with respect to 
the position and the strength of the scissors mode, some of them [7] giving also reasonable 
fragmentation of its strength. However, the situation is more obscure in regard to simple 
phenomenological models whose aim is to explain the physics of the phenomenon and to 
interpret it in the most simple and transparent terms. A noticeable discord of the opinions 
of various authors must be observed here as has been pointed out in [5]. One is forced 
to conclude that there is no general agreement in the understanding of the nature of this 
curious phenomenon. We here will try to shed some light into the confusing situation. 

Our model Hamiltonian will be the one of the well served harmonic oscillator plus 
separable Quadrupole-Quadrupole (Q-Q) residual interaction. The interaction will have 
isoscalar and isovector parts whose coupling constants are reasonably well known from 
the literature. Of course, given such a simple schematic Hamiltonian the RPA equations, 
which are the standard tool to describe the scissors mode, can readily be solved. This, 
however, does not advance us in its physical interpretation. Namely, numerous calcula- 
tions, performed during 25 years of investigation of the scissors mode have undoubtedly 
demonstrated at least one fact: the scissors mode is a very non-trivial, subtle type of mo- 
tion subject to the influence of many factors. In fact its nature is much more complicated 
and interesting than the above picture of two counter rotating and oscillating pieces of 
matter might suggest. First of all one easily imagines that the rotational oscillations of 
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neutron and proton systems are inevitably accompanied by the (quadrupole) distortion 
of their shapes. This entails that the scissors mode is intricately entangled with the Iso- 
Vector Giant Quadrupole Resonance (IVGQR). Thus, if one wants to observe the scissors 
mode, one has to excite at the same time the IVGQR. The quadrupole distortions give 
rise to high lying excitations because of the so-called nuclear elasticity (or Fermi surface 
deformation), the quantum effect discovered by Bertsch [8]. It turns out that without 
this additional restoring force the scissors mode would actually be a zero energy mode! 
Hence, the scissors mode is a pure quantum mechanical phenomenon, which can not be 
explained in the frame of classical mechanics. The above properties of the scissors mode 
highlight perhaps its most characteristic features. They will obtain a natural, sufficiently 
simple, and visual explanation in the frame of our approach outlined below. 

One further issue, strongly debated in the literature, is whether the neutron and proton 
fluids really spatially separate at least to a certain extent during the scissors motion. With 
diffuse surfaces and small amplitude motion this is a not completely trivial question and 
it only makes sense to speak about the separation of the symmetry axes of neutron and 
proton distributions. The debate is not without foundation, since there were attempts [5] 
to construct a model with sharp surfaces: neutron and proton liquids perform out-of-phase 
rotational oscillations with Steinwcdcl- Jensen boundary conditions. Our approach allows 
us to derive the analytic expressions for the current lines and the corresponding figures 
show unambiguously that indeed a separation of the two fluids occurs. More surprisingly, 
the separation not only exists in the low energy magnetic dipole excitation 1+ (the scissors 
mode proper) but also in the K'^ — 1+ branch of the IVGQR, which is named "the high 
energy scissors mode" and whose existence was guessed by various authors many years 
ago [9, 6]. 

From the above discussion it seems clear that in order to elucidate all these subtle 
features an approach involving macroscopic quantities as dynamical variables is indicated. 
Most naturally this can be achieved by working in phase space [10, 11]. This can be 
easily performed by applying the Wigner transform to the Time Dependent Hartree-Fock 
(TDHF) equations. The method of the Wigner Function Moments (WFM) [11, 12] is 
then applied. It can be characterized as a link between microscopic and macroscopic 
approaches: starting from a microscopic Hamiltonian one derives macroscopic dynamic 
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equations for collective variables. Usually one has to establish the set of such variables 
with the help of some physical considerations. The WFM method allows one to avoid 
this non-trivial problem: if one knows at least one collective variable (in our case it is 
the relative angular momentum of neutrons and protons), the procedure of derivation of 
dynamical equations will automatically generate all the other variables needed. 

In this way an unambiguous set of coupled equations in terms of dynamic physical 
variables is obtained, which in the small amplitude limit (RPA) allows for analytic so- 
lutions. Since our equations are written down in the laboratory frame the total angular 
momentum I is, of course, conserved. In this work we study the case without rotation 
and take / = 0. These remarks are important, since microscopic calculations [5] have 
shown that for the results to be reasonable, it is very important to exclude from the wave 
function the spurious component responsible for the rotation of the nucleus, as a whole. 
Such problems do not arise in our approach, because there is no necessity to introduce 
the intrinsic coordinate system. 

The analytical form of our results is very convenient to study the deformation {5) de- 
pendence of various quantities such as position of resonances and transitions probabilities. 
In the small 5 limit we mostly reproduce results already obtained by other authors. How- 
ever, for large 5 we obtain predictions for super- and hyper-deformed nuclei. This area is 
practically not investigated at present. The only investigation within a phenomenological 
model [6] and the only existing microscopic calculation [13] are in rather good agreement 
with our results. In [13] pairing is taken into account whereas we have, so far, not consid- 
ered superfluidity in our approach. However, it is well known that at large deformations 
superfluidity is of little influence on the dynamics and this is also the conclusion in [13]. 
On the other hand at small and moderate 5 the influence of pairing may be appreciable 
and certainly our approach must be generalized to include superfluidity in the future. It 
is for example shown that pairing reduces B(M1) - transition probabilities by important 
factors and that this yields agreement with the experimental deformation dependence [5] . 
In spite of the importance of pairing correlations in nuclei we know from other features 
that it certainly changes results quantitatively but it is too weak in nuclei to yield qual- 
itative changes. For example the moment of inertia usually goes only half way from its 
rigid body value to its irrotational liquid limit (strong pairing case) . We therefore believe 
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that the physical insights we will develop in this paper will stay qualitatively correct, even 
if superfluidity is included in a later stage. As we mentioned already, this shall be the 
subject of future studies. 

The paper is organized as follows. In Section 2 the general outline of the WFM for- 
malism is presented. This formalism is applied to the model of a harmonic oscillator 
potential with Q-Q residual interaction in Section 3: the equations of motion for irre- 
ducible tensors are derived and analyzed, the energies of collective isoscalar and isovector 
excitations are calculated. The method of infinitesimal displacements is used in Section 

4 to find the expressions for the nucleon currents of the different modes and to display 
the respective figures. The magnetic and electric transition probabilities are calculated in 
Section 5 with the help of the linear response theory. Sum rules are analyzed in Section 6. 
The most interesting points in the description and understanding of the scissors mode are 
discussed in section 7. The scissors mode in superdeformed nuclei is considered in section 
8. The summary of the main results and concluding remarks are contained in Section 9. 
To obtain a general impression of the approach and results one can omit the sections 4, 

5 and 6 in a first reading. 



2 Formulation of the method 



The basis of our method is the Time Dependent Hartree-Fock (TDHF) equation for the 
one-body density matrix p'^(ri,r2,i) = {ri\p^{t)\r2) : 



(1) 



where H'^ is the one-body self-consistent Hamiltonian depending implicitly on the density 
matrix and r is an isotopic index. It is convenient to modify equation (1) introducing the 
Wigner transform [14] of the density matrix 



/s s 
d^s exp(-ip • s/h)p^{r + 2 ' ^ ~ 2 ' 

and of the Hamiltonian 



H^{r,p) = / £s ex]){-ip ■ s/h){r + 



(2) 



(3) 



Using (2,3) one arrives [15] at 



^ = - sm - V? • V£ 



where the upper index on the nabla operator stands for the function on which this operator 
acts. If the Hamiltonian is a sum of a kinetic term and a local potential V'^{r), its Wigner 
transform is just the classical version of the same Hamiltonian 



H^^py2m + V^{r). (5) 



Then equation (4) becomes 



Expanding in powers of h leads to 

^ + ^E^'^v,r-Ev.v^^vfr + ^ E v.v,v,ywfv^v^,r - ... = o. (7) 

i=l i=l i,j,k=l 

Now we apply the WFM method to derive a closed set of dynamical equations for dif- 
ferent multipole moments and other integral characteristics of the nucleus. This method 
is described in detail in Ref. [11, 12]. Its idea is based on the virial theorems of Chan- 
drasekhar and Lcbovitz [16]. It is shown in [11, 12], that by integrating equation (7) over 
the phase space {p, r} with the weights Xi^Xi^ . . -Xi^^Pi^^^^ . . ■Pi^_iPi^, where k runs from 
to n, one can obtain a closed finite set of dynamical equations for Cartesian tensors 
of the rank n. Taking linear combinations of these equations one is able to represent 
them through irreducible tensors. However, it is more convenient to derive the dynamical 
equations directly for irreducible tensors using the technique of tensor products [17]. For 
this it is necessary to rewrite the Wigner function equation (7) in terms of cyclic variables 

-±- + -y: i-irp-aVar - E (-i)"v_„y-vsr 

m a=-l 

+|l E (-l)"+'^^'^V_.V_.V_,l--V^,V^,V^,r-... = (8) 

a,u,a=—l 

with 

r+i ^ --j={xi + 1x2) , ro = X3, r-i = --j={xi - 1x2) 

and the analogous definitions for V+i , Vq , V^i , andp+i , po , P-i- The required 
equations shall be obtained by integrating (8) with different tensor products of and Pa- 



3 Model Hamiltonian, Equations of motion, Eigen- 
frequencies 

As outlined in the introduction, the model considered here is a harmonic oscillator mean 
field potential with quadrupole-quadrupole residual interactions. Its microscopic Hamil- 
tonian is 

A 2 1 2 Z N 

i=l fi=—2 i j 

2 Z N 

+ 2'^ E (-1)^{E g2M(r^)g2-M(r,•) + E Wg2-M(r,-)}, (9) 

where the quadrupole operator gf2/i = Y^167r/5r^l2|U and N, Z are the numbers of neutrons 
and protons respectively. The Hamiltonian is of the standard Bohr-Mottelson type [18], 
however, with an interaction coupling protons and neutrons (constant R.) and, separately, 
coupling protons and neutrons among themselves (constant k,). This form is appropriate 
for the description of the scissors mode. The corresponding mean field potentials are 

V^ir,t) = \mu:\' + E {-m^Ql,{t) + -^QUt))<l2-,{^) (10) 

for protons and 

y°(r,t) = ima;V^+ {-lYW2^{t) + RQl^{t))q,_,{v) (11) 

^ fx=-2 

for neutrons. The multipole moments Q2^{t) are defined as 

QUt)= |4p,r}g2^(r)r(r,p,t). (12) 
where /ci{p, r} = 2(27r/i)~^ /ci^p /dV. Introducing the notation 

one can rewrite the mean fields in a more compact way 

V^r,t) = Imu'r' + E {-^rDl,it)q,.,{r). (13) 

^ /i=-2 



/3 ■ 5 

Substituting the spherical functions by the tensor products r^Y2ij, — J — — r^^ , where 



and Ci^i^ is the Clebsch-Gordan coefficient (let us recall that the vector r is a tensor of 
rank one), one has 



Here 

= XR2, + xR% , = xm, + , 

RUt)= I d{p,r}rlj^ir,p,t). (14) 

Integration of the equation (8) with the weights r^^ , {rp)\^ = {r and yields 

the following set of equations (it is important to note, that because (13) is of quadratic 
form, the WFM break off at second order): 



-rRIa Lla = 0, A = 0,2 



dt m 



-^P;^ + 2mc^2L^^-4V^i:y2jTT{^il}(ZJLJ),^ = 0, A = 0, 2 (15) 

3=0 

where {2A1} is the Wigner 6j-symbol. For the sake of simplicity the time dependence of 
tensors is not written out. Further the following notation is introduced 

PUt) = / rK^r(r, p, t), Ll^{t) = / d{p, r}(rp),^r (r, p, t). (16) 

It is necessary to say some words about the physical meaning of the collective variables 
introduced above. By definition K^^ — Q\J -\/6 and is the quadrupole moment of the 
system of particles and i?Qo = — <5oo/v^ Qoo — < > being the mean square 
radius of the same system. By analogy with these variables, defined in the coordinate 
space, we can say that the variables P^^ and Pqq describe the quadrupole moment and 
the mean square radius of the same system in a momentum space. The variables 
describe the coupling of momentum and coordinate spaces. To understand their nature 
it is useful to recall the definitions [11, 12] of nuclear density and mean velocity: 

mn-(r,i)<(r,i) = y J^K/-(r,p,i). (17) 



They enter into the definitions (14,16) of irreducible tensors 

Ll,{t) = jd'rj {rp)x,nv. P,t) = m J dh {ru^h,n^{r, t). (18) 

The last expression for L'^^ demonstrates in an obvious way the physical meaning of 
these variables: being the first order moments of mean velocities they give information 
about the distribution of these velocities in the nucleus. ( "First" means that velocities are 
weighted with the coordinate r). Sometimes, if the motion is comparatively simple, this 
information turns out sufficient to completely determine the velocity field (see section 4). 
In the case of more intricate motions higher order moments are required for a complete 
description of velocities [11]. In any case the moments of velocities are a very convenient 
tool to descibe the collective motion. For example, the zero order moment of velocity is 
nothing more than the linear momentum describing the nucleus' center of mass motion. 
One of the first order moments corresponds to the very well known angular momentum 
of a nucleus. It is connected with the variable by the following relations: 

It is convenient to rewrite the equations (15) in terms of the isoscalar and isovector 
variables 

Rxfj. = Rxn ~ Rxn^ -^A/i = -Pam ~ -^V' -^A/x = -^A/i ~ ^Xfj.- 

So the equations for the neutron and proton systems are transformed into isoscalar and 
isovector ones. The equations for the isoscalar system are given by 

Roo - 2Loo/m = 0, 

Loo - Poo/m + moj^Roo - 2^f5/3[xo{R2R2)oo + Xi(-R2-R2)oo] = 0, 
Poo + 2muj^Loo - 4.^Jb/i[xo{R2L2)oo + Xi(^2^)oo] = 0, 
i?2;u - 2L2^/m = 0, 

- P2jm + muj^R2^ - 2^Jl/i[xo{R2Ro)2^. + Xi{R2Ro)2^ 

-sJV^[X0{R2R2)2^. + Xl{R2R2)2i.] = 0, 
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0)2^ + Xl(-R2-^o)2/i] 

+2V3\xo{R2Li)2^ + Xi(^2Xi)2m] = 0, 
Li, = 0. (19) 

and the ones for the isovector system read: 

^00 - 2Loo/m = 0, 

-^00 - Poo/m + muj^Roo - 2^J^x{R2R2)oo = 0, 

>oo + 2muj''Loo - ^^Jh/i[xo{R2L2)m + Xi(i?2^2)oo] = 0, 

R2n - SZa^/m = 0, 

L2t, - P2jm + muj^R2^ - 2.Jl/i\xo{R2Ro)2^. + Xi{R2Rq)2^] 

-^X{R2R2)2^ - 0, 
4073[Xo(i?2^)2;. + Xl(^2i^ 

-2/773[Xo(i?2l2)2M + Xl(i?2i^2)2^] 

+2V^[Xo(i?2^l)2^ + Xl(^2i^l)2M] = 0, 

^i. + V5x(^2^2)i. = 0. (20) 

Here 

Xo = (X + x)/2 

is an isoscalar strength constant and 

Xi = (X - x)/2 

is the corresponding isovector one. The last equation of (19) demonstrates the conserva- 
tion of the isoscalar angular momentum Li^. The dynamical equation for the isovector 
angular momentum Li^, (the last equation of (20)) describes the relative (out of phase) 
motion of the neutron and proton angular momenta; hence it must be responsible for the 
scissors mode. 

We will need the following algebraic relations: 

(i?2-R2)oO — -^{R22R2-2 — R21R2-I + R20R2O — R2-1R2I + R2-2R22) , 
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2 - - - - 

(-^2-^2)20 = Y ^(-^22-^2-2 + -R2I-R2-1/2 — -^20-^20 + -^2-l-^2l/2 + -R2-2-R22), 

/2 - /3 - 

(i?2-R2)22 = Y ^^(-^22-^20 — V 2-^21-^21 + -^20-^22), 

/2 - /3 - 

(i?2-R2)2-2 = Y ^^(-^2-2-^20 " W 2 "^2-1-^2-1 + R20R2-2), 

[3 - 1 - 1 - 

(i?2-R2)21 = Y y (-^22-^2-1 — ^^-^21-^20 " ^^-^20-^21 + -R2-I-R22), 

/3 - 1 - 1 - 

-{R2-2R2I — -^-^2-1-^20 — -^-^20-^2-1 + R21R2-2): 

(i?2-R2)ll = \j^{.-^R22R2-l — -^-^21-^20 + ^^-^20-^21 " ^^-^2-1-^22), 

(i?2-R2)lO = ^^(y^-^22-R2-2 — -^-^21-^2-1 + -^-^2-1-^21 — Y^-^2-2-R22), 

/s"! - 1- 1-1- 

(i?2-R2)l-l = Y 5 '-^/f"^^^'^^"^ ~ "^-^20-^2-1 + -^-^2-1-^20 " "^-^2-2-^21), 
(i?2-f'l)20 = ^^(-^21-^1-1 — -R2-l-f'll), 

(-^2-^1)22 = ^ -R22L10 — --j=R2iLu, 

/2 - 1 - 

(i?2-f'l)2-2 = Y 2-^2-2-f'lO — ^^-R2-l-f'l-l, 

(-^2-^1)21 = ^^-^22-^1-1 + ^^^-^21-^10 ~ ^^-^20-^11, 
(-^2-^1)2-1 = ~^^-R2-2-^'ll ~ ^^-^2-1-^10 + "^-^20-^1-1- (21) 

With the help of (21) one can write out in detail the whole set of 42 coupled equations 
(including integrals of motion) for the whole set of isoscalar and isovector variables. There 
exists no problem to solve these equations numerically. However, for the time being we 
want to simplify the situation as much as possible what will allow us to get the results in 
analytical form and thus will give us a maximum of insight into the nature of the modes. 

l)Let us consider the problem in the small amplitude approximation, i.e. we only will 
study small deviations of the system from equilibrium. Writing all variables as a sum of 
their equilibrium value plus a small deviation 

Rx,(t) = Rll + n,,(t), P,,(t) = p^,i + V,,(t), L,,(t) = LZ + u,{t), 
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we linearize the equations of motion in TZ\i^, Vx^, and Px^, £'Xn- 

2) We will consider non rotating nuclei, i.e. nuclei with L^^ = L^^ = 0. 

3) Let us consider axially symmetric nuclei, i.e. nuclei with -R2±2 = -^2±i = -^2±2 ~ 

mil = 0. 

4) Finally, we take — Rqq — 0. This means that equilibrium deformation and mean 
square radius of neutrons are supposed to be equal to that of protons. 

After all these simphfications the set of equations for the isoscalar system (19) is 
transformed into the following set of linear equations: 

Ttoo - 2jCoo/m = 0, 

J^oo - Voo/m + ma;^7^oo - 4^1/3x0^20 ^^20 = 0, 
7^00 + 2muj^Coo - ^^XoR2o ^20 = 0, 
ii2±2 - 2C2±2/m = 0, 
il2±i - 2C2±i/m = 0, 
TZ20 - 2jC2o/m = 0, 
J^2±2 - 'P2±2/m + 
Ail - V2±\/m + 
Ao - 7^20/"^ + 



,eq\ 
'20 J 



V2±2 + 2[muj^ - ,j2lixoRll]C2±2 = 0, 

V2±i + 2[ma;' + ^Ji/^x^R2o]^2±i = 0, 

V20 + 2[mw^ + ^xoRVo]^2o - A^XoRVo ^00 = 

jCio = 0. 

/iiii = 0. 



7^2±2 = 0, 
n2±i = 0, 

7^2o-y473xoi?2o^oo = 0, 



(221 



The corresponding set of equations for the isovector system (20) reads 
'^00 - 2jCoo/m = 0, 

>^oo - Poo/m + mu;^noo - ^/^xRll P20 = 0, 



Poo + 2mu;^Coo - 4^1/3xoi?2o -^20 = 0, 
n2±2 - 2C2±2/m = 0, 
p2±i - 2C2±i/m = 0, 
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7^2o - 'iC2o/m = 0, 

^2±2 - V2±2/m + [ma;^ - ^^x^Z - ^Xi^oo 

^20 - V2Q/m + 



^2±2 = 0, 

n2±i = 0, 



ma;2 + ^2/3xi?2^ 



•-00 



4/3xii?i 

^2±2 + 2[ma;2 - .j2/?,xoR7o\£'2±2 = 0, 

^2±i + 2[ma;' + /l/exoi^ao] A±i T VQX0R20 ^i±i = 0, 

ho + 2[mw2 + ^XoRll]C2o - ^XoRll Coo - 0, 



n2o - ^fi/3xon2noo = 0, 



£i±i ± ^j3/2xRlln2±i = 0, 
>Cio = 0. 



(23) 



Due to the approximation 4) the equations for isoscalar and isovector systems are 
decoupled. Further, due to the axial symmetry the angular momentum projection is a 
good quantum number. As a result, every set of equations splits into five independent 
subsets. For example, the equations (23) can be grouped in the following way: 

1) the subset of equations for variables "^oo, >Coo, 'Poo, "^20, ^^20, 'P20 and £10 with pro- 
jections fi = 0. The equation for Cio gives the integral of motion. The rest of equations 
describes the isovector giant monopole resonance plus the branch of IVGQR corresponding 
to fjL — (/9-mode). 

2) the subset of equations for variables "^22, ^^22, ^22 with projections /i — 2 and the 
subset of equations for variables "^2-2, ^^2-2, ^2-2 with projections /i — —2 describe two 
degenerate branches of IVGQR corresponding to /i = |2| (7- mode). 

3a) the subset of equations for variables 'R.21, jC2i,p2i, ^^.u with projections = 1 
describe the scissors mode plus the /i — 1 branch of IVGQR (transverse shear mode). 

3b) the subset of equations for variables "^2-1, ^^2-1, ^2-1, 'Ci-i with projections /i — 
—1 describe the same dynamics as the subset with /i — 1, because the states with fj, — ±1 
are degenerate due to the axial symmetry. 

One also should mention that equations (22,23) are equivalent to the RPA equations 
corresponding to the Hamiltonian (9). The RPA equations have partially been solved in 
[19, 20]. 
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3.1 Isoscalar eigenfrequencies 

The dynamics of the isoscalar angular momentum is trivial - no vibrations, this variable 
is conserved. However it is necessary to treat this mode carefully because, being the 
nonvibrational mode with zero eigenfrequency, it gives nevertheless a nonzero contribution 
to the sum rule (see below) . Let us analyze the isoscalar set of equations with /i — u — 1 
in more detail 



7^21 - 2£2i/m = 0, 



^21 - V2\/m + 



m^2 + ^4/3xoM/\/2-i?^g) 
= 0. 



= 0, 



(24) 



Imposing the time evolution via e'^* for all variables one transforms (24) into a set of 
algebraic equations with the determinant 



in 



-2/m 




-1/m 



2ma;2 + ^2/3x0^20 

in 

The eigenfrequencies are found from the characteristic equation Aj^ = where 



m 



x/3 



(25) 



Using here the relations Xo — 6aco) -R20 — Q2q/V^ and i?oo — —Qqg/V^ (where Qoo — 
A < >) we rewrite it in more conventional notation: 



A,, = 02^02 _ 4^2 _ ^(ge. ^ 4^e.)]_ 



(26) 



For kq we take the self-consistent value (see Appendix): kq — , where uj"^ — -r^,. 



4g, 



00 



Using now the standard definition of the deformation parameters 



we finally obtain 



<52o = Qm\[^^ = Qoo^ ^ 



(27) 
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The nontrivial solution of this equation gives the frequency of the /i — 1 branch of the 
isoscalar GQR 

q2 ^ = 2a;2(l + 5/3). (28) 

In the hmit of small deformation this result coincides with that of [19] . The trivial solution 
= Jlo = is characteristic of nonvibrational mode, corresponding to the obvious integral 
of motion Cu = const responsible for the rotational degree of freedom. Another, not so 
obvious, integral is obtained by a simple combination of the third and first equations of 
(24): 



m 



V21 + 2[muj' + VV6Xoi?2o]ir^2i = const 



V21 + m^a;^(l + -)'R'2i = const. 



Assuming here 6 = we reproduce our result from ref. [12] for spherical nuclei, saying 
that the nuclear density and the Fermi surface oscillate out of phase. 

3.2 Isovector eigenfrequencies 

The information about the scissors mode is contained in the set of isovector equations 
with II — u — 1. Let us analyze it in detail: 



^21 - 2£2i/m = 0, 



jO-21 - V2i/m + 



ma;2 + y^76xi?2g-^/473xli?; 



^21 = 0, 



P21 + 2[mw^ + ^l/6xoi?2o]>^2i - \/6xoi?2o ^11 = 0, 
+ 73/2x^20^21 = 0. 



(29) 



Imposing the time evolution via e*^* for all variables one transforms (29) into a set of 
algebraic equations with the determinant 



in 



>-00 



-2/m 

in 





— 1/m 








in 



/3/2xR2 

Again the eigenfrequencies are found from the characteristic equation Aj^, = where 

A,, = n'[n' - -{muj' + ^Rll)] - -[n'{m^' + ^Rll - %i?^g) - -xxo{R7,f] (30) 
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or in more conventional notation (see the definitions after eq.(25)): 

8 2 Sfi 

A,, ^n^- n2[4a;2 + -n.Qll + -(«:i + 2no)Q2] + ^(«^o - fii)fioiQ2lf- (31) 

m m 

Supposing, as usual, the isovector constant Ki proportional to the isoscalar one, Ki — aKo, 
and taking the self-consistent value for kq we finally obtain 

- 2n^o\2 -a){l + S/3) + - a)S^ = 0. (32) 

The solutions of this equation are 

nl = 0)2(2 - a){l + S/3) ± ^Ju\2 - afil + S/Sf - iu;^{l - a)S''. (33) 

This expression coincides with the result of Hamamoto and Nazarewizh [13] found in 
RPA. The solution 0+ gives the frequency of the /i = 1 branch of the isovector GQR. 
The solution Q_ gives the frequency Qgc of the scissors mode. 

It is worth noticing that in the case Cn — the set of equations (29) becomes quite 
similar to (24). Its characteristic equation reduces to the equation 

- 29w^(2 - a) (1 + 5/3) = 0, (34) 

implying that there exists an integral of motion analogous to the isoscalar one: 

P21 + m^O'^{l + -)^2i = const. 

The nontrivial solution of (34) gives the IVGQR frequency for the case, when rotational 
degrees of freedom are neglected: 

Q2 = 20)2(2 -q;)(1 + 5/3). (35) 

Now let us fix the value of the coefficient a. The experimental fact is: the energy of 
an isovector GQR is practically two times higher than that of an isoscalar one. Assuming 
5 = we have 

Ql = Ql = 2u?{2-a). 

The simple comparison of this expression with (28) shows that the experimental obser- 
vation is satisfied by a = —2. Then equation (33) gives the following formulae for both 
energies: 

El = 4(1 + 5/3 + y'(l7w^)(^^)^ 

El = 4(1 + 5/3 - ^(1 + 5/3)2- ^52 )(^o;)2. (gg) 
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In the limit of small deformations one can write for IVGQR energy 

E,^c.8(l-5/3)(l-^5^)M^. (37) 

For (X——1 formula (35) gives: El^ ~ 8(1 — 5/3)(^t<;o)^. Comparing it with (37) one sees 
that the influence of rotational degrees of freedom on the IVGQR energy is very small. 
The scissors mode energy in the limit of small deformation is 

Esc ^ y^5(l - 6/2)huo ^ ^6huJo, (38) 

which is quite close to the result of Hilton [9]: Esc ~ Vl + 0.66 5 hcuQ. Taking Hcuq — 
45.2/^4^/^ MeV (what corresponds to ro = 1.15 fm used in [21]), one obtains 

Esc ~ 55.4 5A-^/^ MeV, 

which practically coincides with the result of Lipparini and Stringari [21]: Esc — 56 5A~^^^ 
MeV obtained with the help of a microscopic approach based on the evaluation of sum 
rules. Both results are not very far from the experimental [4] value: Esc ~ 66SA~^/^ MeV. 

We now will present the calculation of the flow patterns but the impatient reader may 
jump to section 7 for further discussion on various aspects of the eigenfrequencies. 



4 Flows 

The flow distributions will be calculated with the help of the method of inflnitesimal 
displacements. This method is based on the rules of variation of integral quantities of the 

object [16]. Its detailed description can be found in [11]. 

Small variations TZ]^^ = 5R\^ and = 5L\^ are naturally expressed in terms of 
variations of rf{v,t) and u'l{v,t) (see deflnitions (17, 18)): 

^I,W= Jd'rrl6n^{r,t), 

J^U^) ^mjdh [{rul^)xM'' + irSu-'h^nl^] ^mjd'r (r^ii^AM^,- (39) 

In the last equation we have supposed that u^^ = 0, i.e., there is no motion at equilibrium. 
The variations Sn and Sui are not independent. A relation between them is obtained by 
means of the continuity equation [16] 

i=i ^ 
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where ^j(r, t) = dxi is an infinitesimal displacement. Let us represent it in the form of 
the series: 

3 3 3 

^i{r,t) ^Gi{t) +J2'^iAt)^3+ J2 ^hjk{t)xjXk+ J2 Clijki{t)xjXkXi ^ (40) 

J=l j,k=l j,k,l=l 

For further use we will conserve only the second term of this infinite series, neglecting 
the rest. This procedure is well founded as explained in [11], so we repeat the most 
important arguments very briefly. First of all, it is necessary to notice that due to the 
triplanar symmetry of the equilibrium shape of the nucleus only the tensors Gjj... with 
an even number of indices will survive after integration over the volume. Further, the set 
of dynamic equations (15) for the second rank tensors allows us to describe only rather 
simple types of motion with ~ Xj. To describe a more refined motion with ~ , one 
is forced to consider the dynamic equations for the fourth order moments of the Wigner 
function (the tensors of rank four). There is a one-to-one correspondence: the more the 
motion is complicated, the larger is the number of moments which must be considered. 

So we take ^J{r,t) = J2^=iG'lj{t)xj. It is convenient to introduce the "cyclic" combi- 
nations of analogously to the cyclic variables in (8): 

and to write them as pji{r,t) — J2t=-i{~^y ^Ji,-u{'t)''^i^- Then 

3 +1 Q T +1 

sn-^ = - Ev,(n^en = - E (-i)'^v.(nv-j, K = ^ = E (-i)'^^;;,-.^^.- 

i=l u=-l v=-l 

Using these expressions one finds 

= - / dhY.Cl!;,,.r^r^ i: (-l)^V^(nV-^) 

cr,v <t>=-l 

= 2E E C'^liA-^fSl^-AhMeq). 

k,K (pjCr,!/ 

Now taking into account the axial symmetry {k — 0) one gets 
2 1 

Exactly the same derivation for jC\^ leads to the following result: 
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— ( 1)'^ ^[("^-^20 -^Oo)^1m,10'^/^,0 ( ^-^20 + -^Oo)(^l/i+l,l-l'S'/^+l -1 + '^lM-l,n'^/i-l,l)]- 

We are interested in '^21 and £11. Remembering that Rqq — —Qoo/V^, R20 — (|)^Qoo<^ 
and Qqq — ^Qoo (due to approximation 4)) we find 

^2±i = ^Qoo[(l - ls)So,±i + (1 + p)S±i,o], 

m 2 - 4 - 

-^lii = ^:y|^oo[(l - g'^O'S'o.ii - (1 + 2'^l'S'±i,o], 

where Sa,^ — S^^^ — S^,^ (and Sa,^ — S^^^ + S^^^). Having in mind the e*^* time dependence 
(vibrational motion), we can substitute Sa,i, by iD,Sa,iy- Solving these equations with 
respect to Sa,u, we have 

Now we use the set of equations (29) to find that £11 = — ^ma;^5(l — a)'R.2i and, as a 
result. 

Introducing the notation 

A = -|[i_2g(l-a)5]/[goo(l-^5)], 

S = + 2^(1 - a)5]/[Q,o{l + ^<5)], (41) 



we finally get 

'S'0,1 = ^"7^21, 'S'1,0 = BTZ21- 
A similar analysis of "^2-1 and allows us to write immediately 

'S'0,-1 = ATZ2-1, S-ifi = BTZ2-1- 

So we have for the "cyclic" displacements: 

P+i = -51,0^0 = Bn2iX3, p_i = ^_i,oro = Bn2-iXs, 

Po = -S'o,ir_i - 5*0 -ir+i = V2A{Ji3Xi + ^232:2), 

where J13 = (^2-1 - '^2i)/2, J23 = ^('^2-1 + ^2i)/2. By definition, the variable is 
a small variation of the tensor J^j — Jd{p,r}xiXjf'^{r,p,t). Cartesian displacements are 
found by elementary means: 

- 1 - - i 

6 = ~^^P-'^ ~ P+^'l ^ "^BJi^xs, 5 = ~^^P-'^ P+^'' ^ V2BJ23X3, 
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Let us analyze the picture of displacements in the plane X\ = ^ (the picture in the 
plane 0:2 = must be exactly the same due to axial symmetry) . Knowing the infinitesimal 
displacements 

1-2 = dy^ V2BJ23X3, iz = dz^ ^/2Aj2zX2, (42) 
we can derive the differential equation for the fiow 

-7- = —r— — > ydy -zdz — 0. 

dz Ay ^ ^ A 

Integrating this equation we find 

y^ z^ 

y^ + — const = c — > 1 7— = 1, 

c c/p 

where f3 = —B/A. Depending on the sign of j3 this curve will be either an ellipse or a 
hyperbola. So, it is necessary to study carefully the structure of the coefficient (3. It is 
convenient to study the coefficients A and B separately. 

Let us investigate at first the case of the scissors mode. Taking Q, — and a — —2 
we have 



A = ^ ^ ' \ ^ . (43) 

V2goo(i - i<^)(i + V3 - yo+W^) 

The bounds of possible values of 5 are determined by the natural requirements ij^l^y^z > 0- 
They give (see Appendix): 

-\<^<l- (44) 



It is easy to check that inside of these bounds the square root ^(1 + 5/3)^ — is real 
and the denominator of expression (43) is always positive. The sign of the numerator 
depends on the sign of 5. Elementary analysis of the function 



F{5) = l + 5/?>-^{l + 5/2>Y-\5^-\5 

shows that F{5) = at 5 = and its derivative F' [5) is negative in this point. This 
means that F{5) > for 5 < and F{5) < for 5 > 0. Hence, A > for 5 < and 
^4 < for 5 > 0. Analogous analysis of 5-dependence of B shows that 5 < for 5 < 
and S > for 5 > 0. So, we can conclude that /3 > for any 5 and the currents in the 
case of the scissors mode are described by ellipses. 
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Let us consider the hmit of small S for an illustration. We have 

1-15 1 + 6/3 - ^(1 + 5/3y - + IS 



1 + 1'^^ 1 + d/3 - 7(1 + 5/3)2 _ 3^2 _ IS 



So, for 5 > the short semiaxis of the current eUipse is = c and the long semiaxis is 
Z2 = c//3~c(l + |5). The eccentricity of this ellipse is 



"cur ^ 2 



which must be compared with the eccentricity of ellipsoid corresponding to the shape of 
the mean field: 



'3" 

Thus, the field of currents and the shape of a nucleus are described by prolate ellipsoides, 
their long (and short) semiaxes being disposed along the same coordinate axis. Figure 1 
illustrates the situation schematically. The displacements and the difference in eccen- 
tricities are exaggerated on purpose to demonstrate clearly the essential features of the 
motion corresponding to the scissors mode. One can easily see that its main constituent 
is a rotation (out of phase rotation of neutrons and protons). It is also seen that the 
rotation is accompanied by the distortion of the nucler shape - at least it is evident that 
the long semiaxis becomes smaller. 

Hence, the real motion of the scissors mode is a mixture of rotational and irrotational 
ones. To get a quantitative measure for the contribution of each kind of motion, it is 
sufficient to write the displacement ^ as the respective superposition [5]: 

{ = ae; X f + hV{yz) = a(0, -z, y) + 6(0, z, y). 

Comparing the components = {b — a)z, = {b+a)y with ^2, ^3 i^i (42) we immediately 
find 

b-a^V2J23B, h + a = V2J2zA — > a = D(l + /3), 6 = D(l-/5), 

where D — J2zA/\/2. So in the considered example we have 

a = 2L>(l--5), 6=-L>5, h/a^-8{l + -5)^-5. 

4 2 AAA 
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Figure 1: Schematic picture of isovector displacements for the scissors mode. Thin eUipses 

are the Uncs of currents. The thick oval is the initial position of the nucleus' surface 
(common for protons and neutrons). The dashed oval is the final position of the protons' 
(or neutrons') surface as a result of infinitesimal displacements shown by the arrows. 
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This value of the ratio h/a has the same order of magnitude as another measure of an 
"admixture": the ratio B{Ml)i^/ B{Ml)sc ^ ^5 (see the subsections 7.2 and 7.4). We 
also have to mention the following interesting fact. As we know (see section 7.3), in the 

absence of the Fermi Surface Deformation (FSD) the scissors mode is a zero frequency 

mode. Calculating (with the help of formulae (41)) the ratio B/Aioril — Q we find that 
1 - ^5 

(3 — j^l^ ~ 1 — 2(5. Hence, the eccentricity of the current ellipse (e^„^ = 25) coincides 
with that of the body eUipsoid. As a result, the hues of flows are tangential to the nuclear 
surface, i.e. the motion goes without any separation of neutron and proton surfaces (in 
agreement with the results of papers [22, 23]). Looking at Fig. 1 we can conclude that 
the inclusion of FSD inevitably leads to the separation of neutrons and protons, what 
justifies the name of scissors mode, independently of how large the separation is. 

Let us investigate now the structure of flows for the high- lying mode (IVGQR). 
Taking in (41) Q = Qiy and a = —2 we find 



V2Qoo(l - lS)il + 5/3 + ^(1 + 5/3)2-1(52) ■ 

It is obvious that the denominator is positive in the above-mentioned bounds: — | < (5 < |. 
Elementary calculations show that the numerator is equal to zero at 5 — 3/2, being 
positive for (5 < 3/2 and negative for 5 > 3/2. Hence, A > ior 5 < 3/2 and A < ior 
5 > 3/2. An analogous analysis of the expression 

^ 3(1 + (5/3 + ^(1 + (5/3)2 -1(5^ + 1(5) 

^/2Qoo(l + IS) (1 + 5/3 + ^(1 + (5/3)2-1(52) ' 

shows that B is equal to zero at 5 = —3/4, being positive for 6 > —3/4 and negative for 
5 < —3/4. So we can conclude that < for —3/4 < 5 < 3/2. Hence, the currents in 
the case of IVGQR are described by a hyperbola. As usual, the situation is illustrated 
for the case with small 5. We have 



^ ^-¥ l + S/3 + ^{l + 6/3)^ + i6^ + l6 ,ew,^3., ,^ 1, 

= TT , ^ -(1 - 2d)(l + -5) pa -(1 5). 

1 + 1 + (5/3 + ^(1 + (5/3)2 - |(52 - §(5 ^ ^^2^ ^2^ 

The family of curves y"^ — {1 — \5)z'^ = c is displayed schematically in Fig. 2. The most 
remarkable property of the current lines is seen with one glance: they are nonclosed, 
demonstrating the typical sample (see, for example, first page of [24]) of irrotational 
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Figure 2: Schematic picture of isovector displacements for the high- lying mode (IVGQR). 
The lines of currents are shown by thin curves (hyperbolae) . The thick oval is the initial 
position of the nucleus' surface (common for protons and neutrons). The dashed oval 
is the final position of the protons' (or neutrons') surface as a result of infinitesimal 
displacements shown by the arrows. 
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motion. Nevertheless, the final position of the surface of the proton (or neutron) system 
looks, as if this system was rotated on the whole, the length of the big semiaxis being 
increased. That is, from the outside one again sees practically the same picture, as in the 
case of the scissors mode: rotation plus distortion. This is a curious property of the shear 
motion and it justifies the second name of the — 1"*" branch of IVGQR as "the high 
energy scissors mode" [9, 6]. The quantitative contributions for the two kinds of motion 
to the IVGQR are 

2 ^ 4 ' 4 ^ 4 ' 4 

Concluding the comparison of the scissors mode current with that of IVGQR it is 
worth noticing that two principally different types of infinitesimal displacements result 
approximately in the same change of the nuclear surface position. 

5 Linear response and transition probabilities 

A direct way of calculating the reduced transition probabilities is provided by the theory 
of linear response of a system to a weak external field 

d{t) = 6 exp{-int) + 6^ exp{int). (45) 

A convenient form of the response theory is e.g. given by Lane [25] . The matrix elements 
of the operator O obey the relationship 

I < ipa\d\^o > 1^ = ^ lim (O - 0„)< il:'\d\ip' > exp{-m), (46) 

where ipQ and ipa are the stationary wave functions of unperturbed ground and excited 
states; ip' is the wavefunction of the perturbed ground state, fla — {Ea — Eo)/h are the 
normal frequencies, the bar means averaging over a time interval much larger than 
Q being the frequency of the external field 0{t). 

To use formula (46) in the frame of our method, one must solve two problems [11]: 

(1) to express the matrix element < ijj'\0\ijj' > in terms of collective variables of the 
system, 

(2) to find the solution of the dynamic equations for these variables in the presence of 
the external field. 
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The first problem is solved with the help of the formula for the Wigner transformation 
of a product of two operators [15] 

< > = J d^r J d^r'p{r,r',t)d{r',r) (47) 

= I"'^!^ ^''p (|(^? ■ - ■ ^^)) •"^f - p- 

To deal with the second problem we add the field (45) to the mean field potential (13). 
The equation for the Wigner function (4) is then modified by the term 

Fe,, = ^ Sin ^^{V^ ■ - ■ VrO) {Ow exp{-tnt) + O*^ exp{tnt)) f. (48) 

Proceeding in the same way as before one obtains equations for all collective variables 
needed to calculate < ilj'\0\ip' >. The only new element now is the presence of the term 
Fgxt that makes the equations for the moments inhomogeneous. 

5.1 B(Ml)-factors 

To calculate the magnetic transition probability, it is necessary to excite the system with 
the following external field: 

eh 1 

Ox^' = -^— T— T^(^^^V') • [r x V]. 
mc A + 1 



We are interested in the dipolc operator (A = 1). In the cyclic coordinates it looks like 

Oi.' = -^/Je CI^Jm^.^^, 61 = -01^, = i-iyo,.,. (49) 
Its Wigner transformation is 

(Oi^')w = l^Cli^'la^^P<' = l{rp)ln', 

where 7 = —i \\ — . For its matrix element we have 

' 2mc V 27r 



< i^'\6^A^' >= ^L\^, = ^-{U,, - U,,) = |(£i^, - C,,,). (50) 



Here we have taken into account that L^^, = L^^, = 0. The contribution of Oi^/(t) to the 
equation for the Wigner function is 

Fe.t = 7 {f,. exp-^"* +(-l)'^>-M' exp^"*) 
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with 



Integration of F^/ with the weights r^^ , {rp)xfj_ and yields 



fe,7r 



fc.TT 



^{P, r^F,' = 2^3(2A + 1) ^ C7,^;,v{linn^.(eg). 



fc.TT 



A simple analysis of these expressions shows that the external field modifies only the 
proton part of the set of equations (15) with A = 2: 



dt 
d 

dt 



iQt 



We put here L2^(eg) = 0. The modifications of the respective isoscalar and isovector 
equations are obvious. 

The iJ,' — component of the external field does not disturb a nucleus due to its axial 
symmetry. Let us consider the case of /i' = 1. The set of equations (51) reads 

= 0, 



5^_iexp +5^,1 exp 



rP — . 



K + l\/^P2S K,-^ exp-^* +5^,1 exp^^* 



(52) 



We have used herein the relations i?2o(6?) = -R20/2, Pioi^Q) = ^20/'^ which hold true due 
to approximation 4). 

Now, in accord with formula (50), we have to find the tensors £11 and £11. The tensor 
£11 is found by solving the modified (as in (52)) set of equations (29): 



n2i-2C2i/m = -TyS/Si^^gexp^^*, 



C21 - V2i/m + 



ma;^ + Vl/6xi?2o-\/4/3xii?i 



^21 = 0, 



Ai + y372xi?2o^2i = 0. (53) 
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It is clear that the time dependence of all variables must be exp'^*. The required variable 
is determined by the ratio of two determinants 

A. 



Ai = -r^exp^^ 



where Aj^ is defined by (31) and 

3 



2-' 20 



m 

At equilibrium the set of dynamic equations (15) considerably simplify turning into the 

set of equations of equilibrium. Taking into account one of them 

m V3 v6 



we obtain 

3 

4 m 
Looking at the isoscalar counterpart of the set of equations (53) 



Ar = ^jRQU^cu'' + -(6g2o + 8Qoo) - ^^']- 



7^21 - 2C2i/'m 



- V2\/m + 



mu;2 + ^Xo(i?2g/V2-i?^g) 



7'2i + 2[rm^' + Vl/6xoi?2oK2i 

one easily finds that the isoscalar tensor Cu = 0. 
Writing now the determinant Aj^ as 

we easily can find the hmit (46). For the case where |'0a >= j^'sc > we have 

Applying the standard values of parameters 

Ki = aKo, 4:KoQoo = -muj"^, kqQ 



7V3/8i?^^xp^™, 
0, 

7^378^1 exp'^*, 
(54) 
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S 2 

■-moo 



we arrive at a rather complicated function of the deformation parameter 5 

1 — rv 

I < .c|On|0 > r = -^mu^'QooS'lEl - 2(1 + S/3){hu;mEUEl - El)] 



(55) 



(56) 
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eh 

where hn — - — and E^^ and Ef^ are given by (33). For small values of 5 this expression 
is considerably simplified. Assuming a = —2 one finds the formula 

|<.c|On|U>| //x^-y^ie^ ^ 1 + 5/6 V2167r h ^' 

demonstrating the familiar [5] linear dependence on S. 
For the case {i/ja >= \ipiv > formula (46) gives 

1 — (T 

= -^mu'QooS'[El - 2(1 + 6/3){nwmEUEl - E^)] (57) 
For small values of 6 this expression reduces (for a = —2) to 

I I 111 1/^^ 647r h V^(l + (5/2) V2647r h ^ ' 

Exactly the same results are obtained from the set of equations for the variables 

'^2-1, ^2-1, £-2-1 perturbed by the operator Oi_i. 

Taking into account the relation Qqq ^ ~ - (2^) ' ^^^^^ usually [26] used to 

fix the value of the harmonic oscillator frequency ujq , we obtain the following estimate 

for the transition probability of the scissors mode: 

S(M1) t= 2 I < sc|On|0 > I' = (3/2j'^% 4/3^^2^ ^ QM2A^/H fx^, 

which practically coincides with the result of [21]: B(M1) |= 0.043^*^/^5 /x^, obtained 
with the help of the microscopic approach based on the evaluation of the sum rules. 



5.2 B(E2)-factors 

Electric transition probabilities can be found exactly in the same way as the magnetic 
ones. To calculate the B(E2)-factor it is necessary to excite the system with the external 
field operator 

62,' = er'Y2^, = (3rl„ 61^, = 6^^, = {-1^ O^-,' , (59) 

where (5 — e^J^. Its Wigner transform is identical to (59): {02f^')w — /^'^i//'- The matrix 
element is given by 

< iP'\d2M' >= PK' = ^/^(^v - R2^'). (60) 
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The contribution of 02fj,'{t) to the equation for the Wigner function is 



with 



U,(T 

Integration of F^' with the weights r^^ , {rp)\n and yields 

/ ^^{P, ^Vx^F,, = 0, 

j d{p,r}{rp) 



5(2A+l)i:C^A.:2,'{m}^L(eg), 



/ d{p,r}plF,, = [1 + (-1)^]^5(2A + 1) E^5,2,'{m}^L(eg). 
The external field modifies the set of equations (15) in the following way: 



dt 
'^W - 



+ 



(2x/5C°°,,,i?go(eg) + V7C'/^^;;,Rl^^^,{eq)) exp-^* 
(-1)'^' ,_,,i?^o(eg) + v^Cj,2^^>^,_,,(e5)) exp^^*] . (61) 



The /i' — component of the external field does not disturb a nucleus with axial symmetry. 
Let us consider the case oi /i' — 1 {jj! — —1 gives the same result). The equations (61) 
then read 



L% + ---- -,QZ ^,,-1 exp-^^* exp 



JClt 



(62) 



Now, according to formula (60), we have to find the tensors 7^2i and 7^2i- The tensor 
TZ21 is found by solving the modified (as in (62)) set of equations (29) 



- V2i/m + 



mu;2 + yi76xit:^g-V^xiit:; 



7?.2i - 2C2i/m 

^21 



0, 



3ra+4Q2'o)exp^"*, 



P21 + 2[ma;2 + ^l/6xoi?2o]>^2i - V6Xoi?20 '^n 



0, 



21 — 



4Q^gexp^^*. (63) 
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It is obvious that the time dependence of all variables must be exp'^*. The required 
variable is determined by the ratio of two determinants 



K21 — -T— exp , 



where Aj^ is defined by (31) and 



The tensor TZ2i is found by solving the modified (as in (62)) set of equations (24): 

n2i-2C2i/m = 0, 



- V2i/m + 



mu;' + ^XoiI^l/V2-Rll)\n2i = -|(QoU ^CT exp^^*, 
P2i + 2[mcu2 + 076Xoi?2oK2i = 0, 



Ai = fQ2'oexp*^*. (64) 

Again it is obvious that the time dependence of all variables in these equations must be 
exp*^* and the required variable is determined by the ratio of two determinants 

7^2, = ^exp'-^ 

where Aj^ is defined by (26) and At^ = — A^. 

The hmit (46) is calculated with the help of expression (55) for Aj^, and the analogous 
expression for Ajg: 

In the case \ipa >— |'0sc > we find 

\ < sc\d2,\o > \' = -(3^An{n,,)/mi-nl)2n,,] 

^ ^^3^-^^°° + 4^"°^ + i^'^o(g2^)1/[i^.c(i^L - El)] 

^ ^lL^oo[2(^^'^)'-(l + ^)^'c]/[^.c(£;^-^D]- (65) 

For small 5 (and a — —2) 

I i/S in 12 5 ^0 (5 e^h 5 ^0 J 

< sc O21 > r ~ l^Qm--, tit: ~ l^Qm^- 

' I I ma;ol28V67r °°1-V6 ma;o 128^677 °° 
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In the case l^'a >= li'iv > formula (46) gives 
I < i^;|C>2i|0 > I' = -P^A-^{fli,)/[2ni,{fll-fll)] 

= ^^^Qoo[{l + l)El-2{hu;SmEUEl-El)]. (66) 
For small values of S (and a = — 2) this expression reduces to 

< iv\O2i\0 > p ~ ^Qoo:; r77 ~ ^<5oo- 

' ' ' mcjo 32^271 °° 1 + 5/6 ma;o 32\/27r °° 

In the case iV^a >= iV'is > formula (46) gives 



^[g^llQo^O + 4^2^) + |^'^0(g2^0)V[i^.s(Sl - /.^^^^)] 

e^h 5 „ ^2 nfir.- x\'^^ /[TT i3 



m 167r^°°[^^ + 3^^- " 2(fuD5y]m,r. (67) 
For small values of 5 this expression reduces to 

I < is|02iio > r ^ ^-VgSo(i + m- 

mujQ 16\/27r 

Formula (46) allows one to calculate the matrix element | < 'tpa\0\'>po > P also in the 
case when \%l)a >— |Jlo >, i-e., for the rotational state corresponding to the trivial solution 
of (27): 

I < Q0IO21IO > I' = -/3^A^(Qo)/[2Qo(Jl^ - ^l)] 

= ^[^^"^M + ^Qll) + ^^oiQlimmih'nl - El)] 

= ^^QooSV[hno2{l + 6/3)]. (68) 
m on 

The value of this matrix element is obviously infinite due to the zero value of Qq- However, 
below this expression will be useful to calculate the energy weighted sum rule. 
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6 Sum rules 
6.1 Magnetic case 

The magnetic dipole operator (49) is not hermitian. By definition it is a Hnear combination 
of hermitian operators (components of the angular momentum) 

/V % ^ ^ /V % ^ ^ 

This fact allows one to derive several useful relations: 

[On, [H, = ^([4, [H, Q] + [4, [H, 4]]), 

< O|0ii|i/ >< i/|Oi_i|0 > - ^(1 < z/|4|0 > p + I < z/|4|0 > n 

= -(I < z/ioiiio > 1 < z/ioi_iio > n. 

Using these formulae and the standard sum rule for a hermitian operator [18] 

^(E, - Eo)| < v\U\^ > I' = ^ < 0|[/„ {H, /,]]|0 >, 
one immediately obtains the sum rule for 0\±\. 

Y.{E, - Eo)(| < ^^|On|0 > I' + I < z/|Oi-i|0 > |2) = - < 0|[6n, [E, C>i_i]]|0 > . (69) 

V 

It can also be calculated in a more direct way: 

<o|[6ii,[//,6i_i]]|o>= (70) 

= Y^{Ey - Eo){< 0|6ii|z/ >< z/|Oi_i|0 > + < 0|6i_i|z/ >< u\du\0 >) 

= ^(£;^ - £;o)(< o|6ii|i/ >< o|6I_i|i/ >* + < o|c>i_i|i/ >< o\d\i\u >*). 

V 

Using here the hermitian conjugation properties (49) of the operator Oi^, one reproduces 
formula (69). 

The double commutator is calculated with the help of (9) and (49): 

Z TV 

[Oi^, [E, 0,^\\ = E E Y.{-^rc\Uc\t,2AS^)rU3)^\. (71) 
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Taking into account axial symmetry, one finds the ground state matrix element of (71) 
(in Hartree-Fock approximation) 

< O|[Ol0, [H,di^']]\0 > /fJ,% = ^xX)(-l)''C'2^,2oC2-i/,20-R20-R20 = ^'^</'' -<!!'X(C2^,20-^2o)^- 

It is obvious that this expression is different from zero only for = ±1. Hence, the final 

expression for the right-hand side of (69) is 

9 1 — a 

< 0\[On, [/^,Oi-i]]|0 >= —milfl^l = -^Qoomou'S'i^l = -(1 - a)Eo, (72) 

_ 9 

where, for the sake of convenience, the notation Eq = — — Qoo^'^lAf introduced. The 

47r 

left-hand side of (69) is calculated trivially by multiplying the right-hand side of (56) by 
Esc and adding it to the second line of (57) multiplied by Ei^: 

Ht^t = $^(£;.-£;o)(| <HC'ii|o> I' + l <i/|6i_i|o> |2) 

V 

= 2 (Esc\ < sc\du\0 > |2 + Eiy\ < iv\du\0 > 1^) 

= + Ej^ = (1 - a)So , (73) 



where 



and 



_ [El-2{l + 6/3){nu;r] 



(74) 



[^^,-2(l + ()/3)(/;^)^] ,^ 
^iv^ (^2 _ ^2 ) (1-q;)Eo. (75) 

\ fV SC/ 

So, one sees that the sum rule (69) is fulfilled. 

It is useful to estimate the contribution to the sum rule from the scissors mode in 
the small deformation limit. First of all, with the help of formula (33), one evaluates the 
difference 

El-Elc^2{hon2-a){l + ^-). 
With this the contribution of the scissors mode is calculated quite easily: 

E.e - (1 - a)Eo _ ^2j ^° ^^6) 

Neglecting the 5'^-term and taking hcuo = 41/74^/^MeV (such value was used in the papers 
[20, 27]) we reproduce the result of these papers 

^ 1^ (I) ' ^/^^MeV = 1.4 ^f^6^A^%MeV. (77) 
on \2/ 2 — a 2 — a 

For further discussion see section 7.5. 
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6.2 Electric case 



The sum rule for 02±i can easily be obtained by replacing in formula (70) the operators 
Oi±i by the operators 02±i and using the hermitian conjugation properties (59) of the 
operator 02^: 

Y,{E, - Eo){\ < i^|O2i|0 > I' + I < i^|O2-i|0 > n = - < OIP21, [H, O2-i]]|0 > . (78) 

V 

The double commutator is calculated with the help of (9) and (59): 

p2,, [H, O2,']] = -20/?^- £ E C2';,2,'{l2i}rL(^)- (79) 
Taking into account axial symmetry, one finds the ground state matrix element of (79): 

<0|[O2^,[i/,O2^,]]|0> = -20/32^ E ^2Wi2l}^^0 



A=0,2 



= -2/3^-5^,-.^' \{-\f-j^Rl^ + -^i?2'oj • (80) 
Taking here = 1 we obtain the final expression for the right-hand side of (78) 

< OIP21, [i^,O2-i]]|0 >= -2/32- C-Ql^ + = -e2-;^goo(l + V3). 

m V 3 6 / m 47r 

The left-hand side of (78) is calculated by summing expressions (65), (66), (67) and 
(68) multiplied by the respective energies. It is convenient to calculate the isovector and 
isoscalar contributions separately. The contribution of the isovector modes is 

Q2y^ 1 

2(e,c| <sc|O2i|0> I' + ^i^l <i^;|O2i|0> I') = ^(<3oo + 7(320) 

= e'-^Qoo{l + Sm. (81) 

Exactly the same result is obtained for isoscalar modes: 

4-2 f- 

2 (hflol < 5^0102110 > |2 + Eis\ < is\d2i\0 > \^) = e' — — Qoo(l + S/3). (82) 

Hence the sum rule (78) is fulfilled. 

It is interesting to compare the contributions of the scissors mode and the rotational 
mode. The scissors mode (for small S) yields: 

5 

2Esc\ < sc\O2i\0 > I' ~ ——e^—QooS''. (83) 

1287r m 
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The rotational mode yields: 

2^no| < no|02i|0 > |2 = ^e2-<3oo:r^^. (84) 

Stt m 1 + 0/6 

It is seen that the contribution of the rotational mode is approximately 16 times larger 
than the one of the scissors mode. This is a very significant number demonstrating the 
importance of excluding the spurious state from the theoretical results. Indeed, to describe 
correctly such a subtle phenomenon as the scissors mode, it is compulsory to eliminate 
the errors from spurious motion whose value can be an order of magnitude larger than 
the phenomenon under consideration. 

7 Discussion 

7.1 Hierarchy of variables 

Let us analyze carefully the set of equations (29). It contains a minimal set of variables re- 
quired to describe the discussed phenomenon - scissors mode. The information of the first 
equation is more or less trivial: the tensor is just the time derivative of the quadrupole 
moment Thus, one can say that equations (29) describe the coupled dynamics of the 
angular momentum jC,ii{t), the quadrupole moment "^21 (^) and the quadrupole kinetic 
energy tensor V2i{t). And, what is of principal importance, the angular momentum does 
not play the key role in this ensemble. It is possible to neglect this variable without 
any serious consequence for the rest of equations, which will in such a case describe the 
isovector GQR (see formulae (34,35)). 

The variables 1Z2i{t) and p2i{t) are of considerably more fundamental character. It 
is obvious that one cannot neglect the quadrupole moment 'R-2i{t) because it is the basis 
of the ensemble and the whole problem loses any physical meaning without TZ2i{t). The 
kinetic energy tensor p2i{t) is responsible for the Fermi surface deformation and must 
be taken into account to correctly describe the elastic properties of nuclei [8] and, as a 
result, to get the correct value of the GQR energy. Thus, one arrives at the conclusion: 
it is impossible to construct a reasonable model of the scissors mode with only one pair 
of variables, the angular momentum £11 (t) and its canonically conjugate variable. The 
scissors mode can not exist independently, on its own, without being coupled to the 
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IVGQR (see below). This conclusion is in absolute contradiction with the original idea 
of N. Lo ludice and F. Palumbo [2] and especially with their two rotor model (TRM) 
underlying this idea. 



7.2 Rotation due to vibration 

Rather soon it was understood [21, 28] that the rotational motion must be accompanied 
by the isovector quadrupole vibration (shear mode), the second kind of motion being a 
small admixture to the first one. However, this statement, being true in essence, can be 
misleading to capture the phenomenon. Indeed, one can easily come to the conclusion that 
the rotational motion is the principal constituent of the phenomenon, and the vibrational 
motion is accessory and can be neglected to simplify the description of the problem, if 
one is interested in qualitative results only. It is easy to see, however, that the real 
situation is exactly the inverse! Our analysis of the set of equations (29) has shown that 
the rotational motion (the variable jCii{t)) exists only due to the vibrational one (the 
variables 'R,2i{t), V2i{t)). If one wants to observe the scissors mode, one has to excite 
the IVGQR simultaneously. The IVGQR can exist without the scissors mode, but the 
scissors mode cannot exist without the IVGQR! Neglecting the coupling to the quadrupole 
deformation in the last of equations (29) would induce a full free counter rotational motion 
of neutrons versus protons! 

The only characteristic, in which the rotational motion exceeds the vibrational one, is 

the value of the B(Ml)-factor (see, however, the section 8). The ratio — 7 ~ ^ 5 

is of the same order of magnitude as the coefficient serving to measure the contribution 
of the vibrational motion to the scissors mode: rj — — ~ - (see papers [28, 9]) and 

in paper [21] (^2/) is a frequency of a giant dipole resonance). 

7.3 Fermi surface deformation 

As a matter of fact, the most important ingredient to the scissors mode is the Fermi 
surface deformation. It can be understood by analyzing formulae for the eigenfrequences 
of all three modes. Neglecting the variable 7^21 in (24) we find the following expression 
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for the frequency of ISGQR: 

TTli a/ o TTl o 

- 2 

For self-consistent value of the strength constant kq — — one obtains — 0. This 

4(^00 

result is quite natural, because the pure geometric distortion corresponding to 7?.2i can be 
produced by the proper rotation of the nucleus, without any disturbance of its internal 

structure. Neglecting the variable V2i(t) in (29) we find that the frequency of IVGQR 
(being determined mainly by the neutron-proton interaction) is changed not so drastically: 

Comparing this formula (for a — —2 ) with (37) one sees, that Q?, ~ ScOq changes to 
~ 6uJq. One should recall that also for the Isovector Giant Dipole Resonance the 
distortion of Fermi sphere plays only a minor role. 

It is also easy to see that omitting p2i{t) in (29), one obtains zero energy for the 
scissors mode independent of the strength of the residual interaction. 

Thus, the nuclear elasticity discovered by G.F.Bertsch [8] is the single origin for the 
restoring force of the scissors mode and also for the ISGQR in our simple Hamiltonian of 
a harmonic oscillator with Q-Q residual interaction. So one can conclude that this mode 
is in its essence a pure quantum mechanical phenomenon. This agrees with the conclusion 
of the papers [22, 23]: classically (i.e., without Fermi surface deformation) the scissors 
mode is a zero energy mode. 

7.4 Deformed oscillator and isovector Q-Q interaction 

It is known that the deformed harmonic oscillator Hamiltonian can be obtained in a 
Hartree approximation "by making the assumption that the isoscalar part of the Q-Q 
force builds the one-body container well" [9]. Thus, neglecting the isovector part of the 
Q-Q residual interaction, i.e. assuming ki — ^ k — R — kq, we have to reproduce the 
known results in the deformed harmonic oscillator model. The formulas of other authors 
are obtained, as a rule, in the small 5 limit. For the sake of convenient comparison just 
the same approximation is used here. We have for the energies 
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which coincides with the results of [9] and [29]. The formulae for magnetic transition 
probabilities are 

which coincide with the results of [5]. Possibly they coincide also with that of [29] (as a 
matter of fact, their values are twice larger, but we suppose it is a misprint). For electric 
transition probabilities we find 

^(^2).cTc. — ^Qoo'^, B{E2Uc^ —-—Q,,, - ^, 

in perfect agreement with [5]. The ratios of different characteristics, calculated with and 
without an isovector Q-Q interaction, are 



^sc / 3 ^iv 



V2, 



risc(/tl = 0) V2' fliy{Kl=0) 

B{Ml)sc /3 B{Ml)i, 9 



S(M1),,(ki = 0) V2' 5(M1),,(ki = 0) 2v^' 

B{E2),, _ 1 B{E2)i, _ 1 

B{E2)sc{ki = 0) 2x/6' B{E2)i^{Ki = 0) 
As we can see, the inclusion of an isovector Q-Q interaction increases the energies and 

B(Ml)-factors of the scissors mode and IVGQR and decreases their B(E2) -factors. It is 

also necessary to emphasize the following important result: the ratio 

RM = B{Ml)i^{Ki = 0)/S(Ml),e(«i = 0) 

coincides exactly with the "admixture" coefficient rj introduced by Hilton [28], what sup- 
ports our idea that RM can serve as a measure for an " admixture" . It is seen that taking 
into account the long range correlations (isovector Q-Q forces), one increases RM by 
the factor 3\^/2. And finally, the ratio ilgc/^scif^i = 0) is quite close to the number 
^(1 + 0.66) found in [9]. 



7.5 Sum rule 

Our discussion of the magnetic sum rules will be based on table 4 of the review of Zawischa 
[5] where the results of different models and approaches for ^^^Dy are hsted. For the sake 
of a convenient comparison it is reproduced here (together with its legend) as table 1. 
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Table 1. The energy- weighted orbital Ml sum rule ExB{Mlorh) T (in units of //^MeV) 
in different models evaluated for ^^^Dy as an example, compared with the values obtained 
from the expressions given by Lipparini and Stringari [21]. The total sum rule strength of 
the model and the part exhausted by the low-cncrgy mode arc given. The schematic RPA 
results are from Bes and Broglia [20], for the RPA with Migdal force entry the data of 
Zawischa and Speth [23] are used. In hues 4 and 5 the deformation parameter 5osc — 0.258 
has been taken. 





Total 


Low energy 


TRRM 


140.8 




Classical fluids 


140.8 





NFD model 


141.6 


37.8 


Deformed harmonic oscillator 


31.0 


15.4 


Schematic RPA 


>70 


24.4 


RPA with Migdal force 


108.2 


40.6 


IBM-2 (with E^^?, MeV) 


12.2 


12.2 


Lipparini and Stringari [21] 


145.4 


35.9 



To complete the legend of this table, it should be said that the results of the lines 
1, 2, 3, 4 and 7 are calculated by Zawischa [5], who used hujQ = 46.5/yl^/^MeV (what 
corresponds to Tq = 1.13 fm) and 5 = 0.302 (except the fourth line). Taking the same 
values of parameters one obtains Eq = 42.4//^MeV. In the case of a = — 2 one finds from 
(73) 

^tot = 3Eo = 127.3/x^MeV. 

It is seen from table 1 that this number does not contradict the " Schematic RPA" and is 
in quahtative agreement with the lines "TRRM" (Two Rigid Rotors Model), "Classical 
fluids" , "NFD model" (Nuclear Fluid Dynamics), "Lipparini and Stringari [21]" and "RPA 
with Migdal force" , being exactly in between the last two results. We have to note that 
the result of the small 6 approximation (a; — > ujq, Qoo Qoo) 

Ht^t = 3Eo = 142.6/i^MeV 

is in excellent agreement with "NFD model" and "Lipparini and Stringari [21]" lines whose 
values were obtained in the small approximation also. Such agreement in the last case is 
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especially surprising if one takes into account that the Ml sum rule is model dependent, 
being determined by the neutron-proton interaction, which is quite different in the two 
papers: quadrupolc-quadrupole residual interaction in our case and that of the Skyrme 
type in [21]. This fact confirms that the model Hamiltonian used in this work is realistic 
enough. The exact contribution of the scissors mode to the sum rule (formula (74)) is 

E^c = 30.9/x^MeV. 

The result of the small 6 approximation (formula (77) with hcuo = 46.5/y4^/'^MeV) 

E^c = 35.6/x^MeV 

is rather close to the exact number, being in good agreement with "NFD model" and 
"Lipparini and Stringari [21]" lines. It is worth noting that the ratio Ssc/E^ ~ 1/3 is not 
so far from the value T^sc/^high — 1/4 predicted in [30] on the basis of some theoretical 
analysis of experimental data (their notation "high" means high energy scissors mode). 

In the case of o; = 0, which corresponds to the deformed harmonic oscillator, one has 
(for 5 = 0.258, as in [5]) 

Etot = 31.2At^MeV, E,e = E^^ = IS.e/x^MeV 

reproducing the numbers of the line " Deformed harmonic oscillator" in table 1 . 

Concluding this subsection we have to say that most of the general observations found 

earlier [5] (such as, for example, S^c/Sj^ — 1/3, E^c ~ 5^) are confirmed by our investi- 
gation. 

7.6 Discussion of other approaches to the scissors mode 

After having, as we think, clearly worked out the physics of the low lying scissors mode 
and its interweaving with the IVGQR, one may ask the question about the status of other 
approaches. As already mentioned all other RPA-type approaches [13, 19, 20, 27] have 
from the numerical point of view the same status as our approach (but this concerns 
mostly only the limit of small deformations). Wc also have pointed out that the original 
model of counter rotating rigid rotors [2, 6] does, in our opinion, not at all grasp the 
salient features of the scissors mode. However, the description of the scissors mode also 
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has been attempted with other quite different approaches like IBM (or IBA) , shell model 
calculations, etc. We think that the conclusions concerning IBM (IBA) are very well 
formulated in the review by Zawischa [5] . Therefore, to demonstrate our point of view we 
will simply comment several citations from [5]. 

i) " The original aim of the IBA was the microscopic explanation of vibrations, rotations 
and transitions between the two bosons. Rotational invariance is maintained throughout. 
Starting from a microscopic Hamiltonian, the configuration space is drastically truncated 
by dealing only with nucleons (or holes) in the valence shell, assuming an inert core, and 
approximating correlated pairs of valence nucleons by bosons [31]. (Mostly monopole and 
quadrupole bosons are considered, for special purposes bosons with angular momentum 
different from or 2 are needed too.)... 

The bosons are treated as elementary units, but the internal structure of the fermion 
pairs they represent, reflects itself in the parameters of the model Hamiltonian which 
are fitted to the low energy spectra and transition probabilities [32]. By construction, 
the model is only suitable to describe low-energy states. To describe the giant dipole 
resonance, p and / bosons have been included ([33] and references therein.) The gi- 
ant quadrupole resonances are outside of its scope, as they involve /S.Nosc — 2 excita- 
tions {Nose is the oscillator quantum number). Wc have seen that (in the scmiclassical 
model, in the deformed harmonic oscillator model and in microscopic RPA) a considerable 
amount of strength is in a high-energy mode-the |i^| = 1 component of the isovector giant 
quadrupole resonance. Due to truncation of the configuration space, this strength is miss- 
ing in the IBM-2 sum rule. We have also seen that (in microscopic RPA) a considerable 
part of the Ml strength resides in the region between 4 and 10 MeV, in two-quasiparticle 
type excitations-all these are not included in the model space of the IBA.... Thus, the 
IBA sum rule by the basic assumptions of the model comprises only a small part of the 
full sum rule. Therefore, it is not so well suited to assess the collectivity of a state." 

ii) "As already mentioned, among the predictions, the IBA was closest to the strengths 
to be detected. The reason is, of course, that the model parameters arc extracted from 
low-energy data and the model is best adapted just to the energy range up to a few MeV 
where the Ml states have been found... The naive assumption of bare orbital g-factors 
already gives a good prediction of the low-energy strength. Adjustment of the model 
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parameters to the low-energy vibrational and rotational states [32] finally has the effect 
that the average energy is also well reproduced." 

iii) "In order to find out the physical interpretation of the lowest K = 1 mixed sym- 
metry state of the IBM-2, the (semi) classical limit of the IBM-2 Hamiltonian has been 
investigated by Balantekin and Barrett [34], Bijker [35] and Walet [36], yielding a Hamil- 
tonian similar to that of TRRM, the potential being a function of the angle 9 between 
the symmetry axes of protons and neutrons: V{9) — X9'^." 

With respect to iii) we must note, that deriving the classical limit, the vibrational 
degrees of freedom have been neglected to simplify the derivation. Not surprisingly they 
obtained the TRRM Hamiltonian as a result for which we have already given our opinion 
above. 

Generally speaking these above citations give exaustive characteristic of the status of 
IBM (IBA) calculations of the scissors mode: they are able to give correct values for the 
energy and strength of the scissors mode but they do not explain the real physics of the 
phenomenon. 

The situation with shell model calculations is rather comphcated. There are the 
well known difficulties to treat heavy nuclei because of the huge dimension of matrices. 
Therefore the calculations are usually made for very light nuclei. Even there one must 
divide them into two groupes. There are qualitative estimations with truncated basis 
{ANosc = 0, see for example [37, 38]). Naturally, those calculations suffer from the same 
drawbacks as the IBM calculations and then the same comment given above apply. There 
are also 'reahstic' calculations (for ^Be and ^°Be) with an extended {ANosc — ^)+{^Nosc — 
2) basis [39]. In principle they have the same status as RPA calculations. Still one can 
ask the question whether it makes sense to talk about scissors mode in such light nuclei. 

We also would like to comment on the neutron-proton deformation (NPD) model. 
We again cite [5]: "The Bohr-Mottelson model has been generalized to isovector degrees 
of freedom leading to the neutron-proton deformation (NPD) model... Rohozinski and 
Greiner [40] applied the NPD model to the orbital magnetic dipole excitations.... 

The TRRM also has been changed by its authors to the TRM, relaxing the condi- 
tion of rigid rotation and replacing the rigid-body moment of inertia by a smaller one 
obtained from some model or a phenomenological one [41]." If their parameters "...are 
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adjusted to the data, the TRM and NPD model will coincide, even though in their orig- 
inal assumptions-rigid rotation in the one, irrotational flow in the other-the models are 
contradictory." 

Therefore wc agree with the conclusion of Zawischa [5] that the NPD model can not 
describe the scissors mode. The reasons are also especially well born out in terms of 
our method: the collective variables in their model are only "^2// and C,2jj. , there are 
no rotational degrees of freedom {Cii,) and there is no Fermi surface deformation {V2n)- 
As we know from an earlier discussion, only the isovector giant resonance survives in 
such conditions (see section 7.3). Therefore, the (5^-dependence of i?(Ml)-factor which 
the NPD model predicts and which is in principle in agreement with experiment, must 
actually be interpreted as the 5^-dependence of B{M1) for the IVGQR (see eq. (58)). The 
5-dependence of the low lying scissors becomes quadratic only after inclusion of pairing 
[5]. 

In conclusion we think that above set of citations and argumentations is convincing 
enough to state that all the models and methodes describing the scissors mode without 
couphng to IVGQR are pure phenomenological and are therefore of restricted usefulness. 

8 Superdeformation 

As already mentioned, a certain drawback of our approach is that, so far, we have not 
included superfluidity into our description. Nevertheless, our formulas (36, 56, 57) can 
be successfully used for the description of the superdeformed nuclei, where the pairing 
is very weak [6]. For example, applying them to the superdeformed nucleus -"^^^Dy (5 ~ 
0.6, hujo = Al/A^/'^MeV), we get 

Ei^ = 23.6 MeV, B{Ml)i^ = 15.9 ix% 

for the isovector GQR and 

Esc = 5.4 MeV, B{Ml)sc = 20.0 i^l 

for the scissors mode. There are not so many results of other calculations to compare 
with. As a matter of fact, there are only two papers considering this problem. 
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The phenomenological TRM model predicts [6] : 

£;,^~26MeV, ^(Ml)^^ ~ 26//^, E^.c^QAMeV, B(M1),, ~ 22//^. 

The only existing microscopic calculation [13] in the frame of QRPA with separable forces 
gives slightly more information: 

- 5 - 6 MeV, B{Ml)i+ ~ 23 iJ,%, B{M1),^ ~ 0.4 

Here B{Ml)i+ denotes the total Ml orbital strength carried by the calculated K'^ — 1+ 
QRPA excitations modes in the energy region below 20 MeV. The B{Ml)sc denotes '"the 
calculated overlap probabilities of the QRPA solutions with the synthetic orbital scissors 
mode which is defined as 

\R>= N-'{CU-QJ^U)\9-s.>, 

where is a normalisation factor. The parameter q is determined by the requirement 
that the mode \R > is orthogonal to the spurious state IS" >~ Ci-i\g.s. >" [13]. 

It is easy to see that in the case of IVGQR one can speak, at least, about qualita- 
tive agreement. Our results for Egc and B{Ml)sc are in good agreement with that of 
phenomenological model and with Egc and B{Ml)i+ of Hamamoto and Nazarewicz. The 
very small value of B{Ml)syn is explained quite naturally by the fact, that the synthetic 
mode does not treat properly two main ingredients of the scissors mode: Fermi surface 
deformation and couphng with IVGQR. 

Examples of ^-dependences of energies and B(Ml)-factors are shown in Fig.3. 

9 Conclusion 

In this work we again have considered the issue of the physics behind the nuclear scissors 
mode. In spite of 25 years of research and many valuable contributions to this subject 
the subtleties of the scissors mode are still under debate, and in our opinion erroneous 

interpretations of the subject continue to appear in the literature. Surprisingly, no sys- 
tematic study of the mode in the Bohr-Mottelson Hamiltonian has yet been carried out. 
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and our purpose here was to fill that gap. The Bohr-Mottelson Hamiltonian consist of a 
harmonic single-particle potential together with a separable Q-Q interaction. The Q-Q 
forces have different couplings in the isoscalar and isovector channels. The isoscalar cou- 
pling strength is determined from Bohr and Mottelson self-consistency condition, leaving 
the isovector strength as a free parameter. We adjust it from the fact that the isovector 
giant quadrupole resonance is experimentally known to lie practically at twice the energy 
of the isoscalar giant quadrupole mode. With this our model is entirely fixed and its 
solution in the small amplitude limit can be found analytically for excitation energies and 
transition amplitudes. 

The physics becomes particularly transparent once the TDHF equations are written 
down in phase space and the so called Wigner function moments are introduced. This 
approach allows one to establish the optimum set of macroscopic variables: quadrupole 
moment, angular momentum, pressure tensor, etc. These variables are, in the scheme 
of our formalism, absolutely unambiguous and, together with the analytic solution, they 
allow for a maximum of physical insight. At least, the inevitable coupling of the scissors 
mode with the isovector giant quadrupole resonance becomes obvious immediately, already 
at the stage of the formulation of the model. Furthermore, the Fermi surface deformation, 
whose decisive role in the physics of the scissors mode is difficult to predict employing 
naive phenomenological models, appears in our approach quite naturally. 

The eigenvalue equation in the isovector channel yields two frequencies which are given 



by n± = 2u^{l + S/3) ± ^(1 + 5/3)2 _ |52_ They are distributed in a non symmetric 
way around twice the harmonic oscillator frequency Cu and they gradually approach one 
another with increasing 6. The low lying frequency corresponds to the one of the scissors 
mode proper, whereas the other is the so-called high lying scissors mode. Our analysis 
shows that, indeed, the motion of both modes is "scissors" -like in the sense that the long 
symmetry axes of proton and neutron distributions get tilted by a small angle during their 
oscillatory motion. Nevertheless both modes are quite distinct what is revealed by looking 
at the respective flow patterns (Figs. 1,2). The flow lines of the scissors mode are closed 
ellipses (i.e. mostly rotational flow) leading to a compression of the long axis, while the 
ones of the high lying mode are open hyperbolas (i.e. mostly irrotational flow) leading to 
an elongation of the long axis. The frequency of the scissors mode as a function of mass 
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number turns out to be about 15% too low, compared with the experimental data. In 
respect to the somewhat crude model we have been employing this may seem a reasonable 
agreement. One should, however, mention that we completely disregarded pairing in our 
work. It is generally known that superfluidity makes the mass parameters smaller, i.e. 
the frequencies of the modes become higher. It will be a further task to study whether 
this accounts for the missing 15% in energy. 

Other quantities we studied in our model are transition probabilities, for instance 
with respect to their deformation dependence. Though for small deformation most of our 
results have already been found by other authors [5] with different methods, we make a 
point here in predicting the behaviour for superdeformed nuclei. 

We also want to attract the attention to the potential richness of the set of our equa- 
tions (15). In a further study one may employ them for the description of the joint 
dynamics of the isoscalar and isovector giant monopole and quadrupole resonances plus 
the scissors mode in deformed rotating nuclei, the amplitudes of vibrations being not 
necessarily small. A large amplitude motion was already treated in the frame of this 
approach to describe the multiphonon giant quadrupole and monopole resonances [12]. 
What about two-phonon scissors? The question is not only academic - the first attempt 
to interpret some numerical results as the multiphonon scissors is already known [42]. 

One may also think to take into account the spin degrees of freedom - only the number 
of dynamic equations must be doubled (spin projections up and down). Then, the theory 
becomes capable of describing spin-fiip excitations. As a result, there appears a possibility 
of considering the orbital and spin components of the scissors mode simultaneously. 

It is worth noticing that the set of equations (15) is written in the laboratory coordinate 
system. It allows one to get rid of any troubles connected with spurious rotation, because 
the total (i.e., isoscalar) angular momentum is conserved (the last equation of (19)). 
Moreover, the total angular momentum enters into the dynamic equations - hence one can 
study the behaviour of all modes in rotating nuclei. This would be especially interesting 
for the scissors mode in superdeformed (SD) nuclei, because "The SD bands in nuclei 
around ^^^Dy and ^^^Hg are observed at high spins" [13]. The first interesting results of 
calculation interpreted as the rotational band built on the scissors excitation appeared in 
[42]. The above mentioned problems shall be investigated in future work. 
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Appendix 

It is known that the deformed harmonic oscillator Hamiltonian can be obtained in a 
Hartree approximation "by making the assumption that the isoscalar part of the Q-Q 
force builds the one-body container well" [9]. In our case it is obtained quite easy by 
summing formulas (10) and (11): 

V{v, t) = hv^ir, t) + V^{r, t)) = \mujh^ + «o E {-lYQ2,{t)q2-,{v). (85) 

In the state of equilibrium (i.e. in the absence of external field) Q2±i = Q2±2 = 0. Using 
the definition [18] Q20 = Qoo^^ and the formula ^20 = 2z'^ —x"^ — y'^ we obtain the potential 
of the anisotropic harmonic oscillator 

with oscillator frequencies 
8Q00 

where a — —kq -. The definition of deformation parameter 5 must be reproduced by 

the harmonic oscillator wave functions, that allows one to fix the value of a. We have: 

Qoo = -(^ + ^ + ^), Q20-2^{^-^), 

where S^. = llf^^{n,j. + -)j and is the oscillator quantum number. Using the self 
consistency condition 

where Eg and ljq are defined for spherical nucleus, we get 

Qoo ^' + 2u;' 3 ■ 
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Solving last relation with respect of a we find 



Now we can write final expressions for oscillator frequences 

2 2 2 2 2 

and the self consistent value of the strength constant 

mcj^ 1 

The condition for volume conservation — const — cUq makes cu 5-dependent: 

2 2 1 + 1^ 

CU — UJn 



°(1 + |5)2/3(1- |5)l/3- 

Qoo depends on 5 as 

Qoo - QSo^^ ^ 4^)1/3(1 _ 2^^2/3' 

where g°o = ^|^^ ^ = ^o^i^/^- 
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